Power fluctuations in the wavelet space: 
large-scale CMB non-gaussian statistics 
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Abstract 



We analyse the large-scale coherence of the CMB anisotropy field with non-gaussian 
initial conditions using 2-point function of the power fluctuations in the wavelet 
space. 

Employing the multivariate Edgeworth expansion (MEE) we constrain the nor- 
\q • malization for the cosmic string mass-per-unit-length /j,, obtaining Gfi/c 2 = 1. 075^0*375 x 

00 ! 1CT 6 at 68% CL in the standard gaussian statistics. This value is consistent with 

the results obtained from the simulations of the evolution of the string network and 
Q ' by other large scale studies. 

00 ' 
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! 1 Introduction 



The cosmic microwave background (CMB) provides one of the most useful 
tools for understanding the physical processes in the early universe, giving 
valuable information about the origin and the dynamics of the primordial 
fluctuations that seeded the today large-scale structures. There are two com- 
peting families of cosmological theories quite distinct in this regards. 
In the models invoking most types of inflation [1-3] the perturbations' gener- 
ation mechanism is linear and gaussian fluctuations are predicted. For a given 
type of inhomogeneity (adiabatic or isocurvature) the only degree of freedom 
in the initial conditions is the amplitude, and since the evolution is linear, so- 
lutions scale linearly with the initial amplitude. Linearity combined with the 
assumed statistical homogeneity of the universe guarantees that one may de- 
compose the cosmological perturbations into eigenmodes with eigenfunctions 
that evolve independently. In the Fourier space, for instance, the evolution of 
a perturbation in each mode k is independent. For a finite sample of the CMB 
field, the observed fluctuations are the convolution of the "true" fluctuations 
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of the infinite field with the Fourier transform of a mask [13]. The 2-point 
function of the the power fluctuations depends only on the selection function 
[13,14]. 

The second class of theories, most topological defect models [4-6] and some 
versions of inflation [7], predict non-gaussian statistics of the primordial per- 
turbations. Those models involve spontaneous symmetry breaking as the hot 
early universe cooled leading to a non-linear perturbations' generation mecha- 
nism. For this type of perturbations the Fourier modes k are coupled and their 
evolution cannot be anticipated on the basis of the initial conditions [8,10]. 
The predicted angular power spectrum has a position dependence with a vari- 
ance larger than the variance of the angular power spectrum for the gaussian 
case at the same multipole order [9,10]. The defect theories are highly con- 
strained by causality, the fluctuations being completely uncorrelated beyond 
the horizon scale at all times [8,11]. 

It is difficult to make difference between gaussian and non-gaussian fluctua- 
tions of the density field at large angular scales because of the tendency of 
any distribution to approach gaussian statistics when averaging over large 
area (the central limit theorem), and also because of the cosmic and sampling 
variances. 

In this paper we propose the wavelet transform to analyse the large scale co- 
herence of the CMB anisotropy field with non-gaussian initial conditions. The 
main advantage of the wavelet transform over the Fourier transform is its ca- 
pability to investigate features of the CMB anisotropy maps with a resolution 
according to their scales [15]. The wavelet base functions are highly localized 
in space and the information regarding the scale and the position are stored 
explicitly. Also, the wavelet analysis can be performed with a finite sample 
of data in contrast with the Fourier analysis where the result depends on the 
selection function. 

We employ 2-point function of power fluctuations in the wavelet space to com- 
pare COBE- DMR 4- year temperature anisotropy maps to the predictions of an 
analytical cosmic string model [21]. The advantage of the analytical approach 
is that one can obtain statistical fluctuations with non-gaussian random phases 
for a given experimental configuration, while the explicit dependences on the 
normalization of the primordial power spectrum, string parameters and angu- 
lar scales can be shown. 

We use the non-gaussian likelihood function to constrain the dimensionless 
cosmic string parameter Gfi/c 2 (where G is the Newton's constant, c is the 
speed of light and fi the the mass per unit length of the cosmic strings) at 
68% CL in standard gaussian statistics, by comparing zero lag autocorrelation 
function from COBE- DMR with the predictions of the cosmic string model. 
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2 Monte-Carlo simulations 



According to the scaling solution for cosmic strings [17,18] there is a fixed 
number M of strings with the curvature radius of the order of the horizon per 
Hubble volume at any given time, with orientations and velocities uncorrelated 
over distances larger than the horizon. The temperature anisotropy induced 
by long strings is given by Kaiser-Stebbins formula [19]: 

ST 

— = ±4nGfi\k ■ (>y s v s x e,)|, 

where: k is the direction of observation, v s is the velocity of the string with 
the orientation e s , 7 S = (1 — Vs/c) 1 ^ 2 , G is the Newton's constant and fx is the 
mass per unit length of the string. 

The time from the last scattering surface until now is divided into N Hubble 
time steps U with t i+ i = 2ti and iV = log2(t /ti s ). The apparent angular size 
of a Hubble volume at time t; is 9 Hi ~ ~ t}^ 3 and 9 Hi+1 = 2 1 / 3 ff . for 

large redshifts in the matter dominated era assuming Q Q = 1. 
For an experiment measuring a square map of size 8° x 6° at each Hubble 
time ti the number of rii string segments random placed over the area size 
(9° + 9 Hi ) 2 is given by [20]: 

n i = M(6P + 6 Bi ) 2 /& 2 Hi . (1) 



We use an analytical statistical string model [21] to take into account the 
combined effects of the temperature fluctuations induced by strings present 
between the redshift z ~ 250 and today (N = 12). 

We calculate the probability distribution functions of the temperature fluc- 
tuations induced by cosmic strings for different values of the string scaling 
solution parameter M by Fourier transforming the characteristic function of 
this model normalized to approach a gaussian distribution with zero mean and 
a = 1 when M -»• oo [21]. 

We generated CMB independent realizations, each of 6144 pixels, using the 
spherical harmonic representation: 

w 30 I 

-T^(M) = E E a /m WW m (M). (2) 

1 1=2 m=0 



where W\ is the DMR gaussian window function with FWHM=7°. 

The harmonic coefficients ai m are random variables with zero mean and the 

variances given by [22]: 

. 2 ^ =n 2 47r r[Z + (w-l)/2]r[(9- W )/2] 

lm ^rms-ps 5 r[/ + (5 _ n )/ 2 ]r[(3 + n)/2] ' 1 ' 
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where: Q rms -ps is the quadrupole normalized amplitude and n is the spectral 
index of the primordial power spectrum. 

The coefficients a\ m are drawn from parent populations obtained in the cosmic 
string model The non-gaussian amplitude distributions of a\ m with random 
phases convoluted with the spherical harmonics Yi m (6, <fi) are characterized by 
a positive kurtosis in the resulting temperature distribution. 
We generated 800 n — 1 full-sky realizations for each model defined by 
Qrms-ps an d M. To each realization we added a realization of the noise deter- 
mined by the instrument sensitivity and the number of observations per pixel 
of DMR 4- year 53GHz (A+B)/2 map, and rejected the pixels with galactic 
latitude -19° < b < 20° obtaining 2 12 pixels. 

We generated also in the same conditions 800 n = 1 full sky realizations for the 
same values of Q rms -ps with ai m drawn from a standard gaussian distribution 
with zero mean and a — 1. 



3 Two point function of power fluctuations in the wavelet space 

A wavelet [15] ty(x) G R is a function whose binary dilatation and dyadic 
translations generate an orthonormal base so that any function f(x) G R can 
be approximated up to an arbitrarily small precision by: 

/(xl^^cj^-j). (4) 

i£N j=0 

The base functions W (x) are labeled by the scale % and the position j G N) . 
The wavelet coefficients c* are real quantities containing all the information as 
that contained by the Fourier coefficients. The major difference is the number 
of indices, that is, in wavelet analysis the scale % and the position j in the 
scale are explicitly stored. The wavelet base functions ty(x) satisfy the scaling 
equation: 

¥ j (x) = 2- i / 2 *(2- i (j-x)). (5) 

Let be the Hubble volume spanned by the scaling function [Q{2 l x — j)]jeN 
at the scale i. From the scaling equation it follows that: 

••• C Vh.^ C V h . C V m+ i C ... 

This hierarchical structure ensures the multi-resolution capability of the wavelet 
analysis. 
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The orthogonality condition: 



( 




) 




(6) 



yields to: 




^f(x k )Hf)(x k ). 



(7) 



k 



This last equation represents our wavelet transform convention. 

For a random gaussian field the ensemble-averaged power P(A6) in the wavelet 

space at the separation angle A9 is: 



< 4(711)4' (rf 2 ) > Ae = 8 i i '5 j j 'P(A9), (8) 



where n[ and ri2 are two directions in the sky separated by A9. 

For a random non-gaussian field generated using the cosmic string model the 

intrinsic correlations exist only within the same scale and the equation (8) 

becomes: 



and Pjj*(A8) contains the correlated ensemble-averaged signal. 

We define the 2-point function of the power fluctuations in the wavelet space 

as: 



In this paper we present the results obtained using Daubechies wavelet base 
functions of order 4 as a typical compact orthonormal base. 
The wavelet transform applied on 2D map (with 2 12 pixels and i = 12 uncor- 
rected scales) acts in a pyramidal algorithm, stripping the initial map scale 
by scale into its correlated components. The reader may refer to the references 
[16,23] for complete information. 

Figure 1 presents the distributions of spots as a function of threshold (ex- 
pressed in number of standard deviations) obtained for the averaged simulated 
maps in real and wavelet space for two types of simulations: a/ m with random 
gaussian phases (panels a and b ), and a\ m with random non-gaussian phases, 
obtained in the cosmic string model with M = 10 (panels c and d ). For all 
simulations we take Q rm s-ps — 18/iK. In each case the gaussian fits are also 
shown. One can see that the differences between the distributions obtained in 
gaussian and non-gaussian case are more evident in the wavelet space than in 




(9) 



(10) 
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the real space. 

The maps obtained by Monte-Carlo simulations were wavelet transformed and 
then 2-point functions of the power fluctuations were constructed (binned in 
n = 36 bins, in the manner of COBE-DMR). 

Figure 2 presents the 2-point function of the power fluctuations in the wavelet 
space £p(A9) normalized at £p(0). The solid curve represents the theoretical 
prediction if the underlying distribution is gaussian and Qrms-ps — 18/xK. 
The error bars are based on Monte-Carlo simulations in the cosmic string 
model with M = 10 and the same normalization. 



4 Non-Gaussian Likelihood analysis 



According to the standard x 2 method, if a specific model defined by the param- 
eters (Qrms-ps, M) is more likely to have occurred higher is the probability 
to obtain values of the x 2 larger than the measured Xo- The cumulative prob- 
ability used to define the confidence regions on the parameters is given by the 
integral: 

Xo 

o 



where L is the likelihood function full specified by the x 2 value and the co- 
variance matrix. We define the x 2 as: 

x 2 = f; f;« «p > -& dms )a-'« a > -6 DM \ (i2) 



where < £p > is the ensemble-averaged value of the 2-point function of the 
power fluctuations for Monte-Carlo realizations and C,p MR is the two point 
function of the power fluctuations for DMR data. 

The covariance matrix \j calculated for the the Monte-Carlo realizations is: 

i N r ealiz 

^ = E < & >M~ < & >)• (13) 



The standard Gaussian likelihood statistic L g associated to the two point 
function of the power fluctuation £ P (A0) is given by: 



e 2 



!„2 



Lg (2vr)™/2( c ietA) 1 /2- ^ 
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In the case of non-gaussian random fields the confidence regions of the param- 
eters depend also on the higher order correlation functions. It is shown that 
for mild non-gaussianity the likelihood function can be obtained within the 
multivariate Edgeworth expansion (MEE) [12]. According to MEE the total 
likelihood function is: 

L — L g + L ng , (15) 

where L ng is the non-gaussian correction which embodies the higher order 
moments of the data. 

The integral in equation (11) can be written as: 

I(xl) = F n (x 2 ) + Q n (xl), (16) 

where F n (xl) is the integral of gaussian likelihood, that is the cumulative in- 
tegral of the x 2 probability distribution with n = 36 degrees of freedom, and 
Q n (Xo) i s a non-gaussian correction that depends on the higher order mo- 
ments. 

For each set of Monte-Carlo simulations we estimated the cumulative proba- 
bility I{ X 2 ). 

Figure 3 presents I(x 2 ) obtained for the models with Q rms -ps = 18/iK and 
few values of the M parameter compared with the same probability obtained 
for the gaussian case (the upper curve). 

One can see that for mild non-gaussian fields the confidence regions are larger 
than the confidence regions in the gaussian case at the same x 2 threshold. 
Using the gaussian likelihood given by equation (14) we found Q rms -ps — 
18.29 ± 1.36/iiT when the DMR 53GHz (A+B)/2 map was compared with the 
Monte-Carlo simulations obtained for a gaussian underlying distribution. 
Figure 4 presents the confidence regions obtained for the parameters Q rms -ps 
and M in cosmic string models at Xo = 31.1 that corresponds to 68% CL for 
a purely gaussian statistics. 

Considering the region with I(xo) > 0.3 and taking the extrema of this contour 
we found that DMR data favourize the cosmic string models with M within 
10 and 12 and Q rms -ps within 17.7QfiK and 20.24/xK at 68% CL for the stan- 
dard gaussian statistics. We estimated the dimensionless cosmic parameter 
Gfi/c 2 by comparing, in the wavelet space, the averaged zero lag autocorrela- 
tion function obtained in those models with zero lag autocorrelation function 
of DMR 4-year 53GHz (A + B)/2 map. We found: 

G/j/c 2 = 1.075±g;p| x 1(T 6 . (17) 

There are two terms that contribute to the error: a symmetrical error (±0.185) 
originating from the width of the confidence region and an asymmetric error 
given by the cosmic and sampling variances. 
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5 Conclusions 



In this paper we propose the wavelet analysis for studying the large-scale CMB 
correlations induced by the cosmic strings. 

Using an analytical cosmic string model to generate CMB temperature fluctua- 
tions with non-gaussian random phases in the COBE- DMR experimental con- 
figuration, and 2-point function of the power fluctuations in the wavelet space, 
we constrain the normalization of the primordial power spectrum Q rms -ps and 
the parameter of the cosmic string scaling solution M. 

We found that the large-scale CMB anisotropy favourizes the cosmic string 
models with large number of strings per Hubble volume, M=10, 11 and 12 
and Qrms-ps from 17.76/iK to 20.24/iK at 68% CL for the standard gaussian 
statistics. 

By comparing zero lag autocorrelation function of COBE- DMR 53GHz (A+B)/2 
map with our predictions we obtained the dimensionless cosmic string param- 
eter Gfi/c 2 = 1.075^^375 x 10 ~ 6 , value that is compatible with the values 
obtained from the simulations of the evolution of the string network from 
z = 100 to present [24], as well as with other existing studies of the expected 
large-scale CMB anisotropy [21], [25,26]. 
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Figure Captions 



Figure 1. 

Number of spots as a function of threshold for an averaged map obtained for: 
a) aim with gaussian random phases, b) the wavelet transformed of the gaussian 
maps, c) a/ m with non-gaussian random phases (see text), d) the wavelet trans- 
formed of non-gaussian maps. The gaussian fits for each case are also shown. 



Figure 2. 

The 2-point function of the power fluctuations in the wavelet space: with the 
continous line is the theoretical prediction if the underlying distribution is gaus- 
sian and Qrms-PS = 18//K. The error bars are based on Monte-Carlo simula- 
tions in the cosmic string model with M = 10 and the same normalization. 



Figure 3. 

I(x 2 ) f° r the models with Q rm s-PS = 18//K and M=4, 6, 9 and 10 (from bot- 
tom to top) and the same probability in the gaussian case (the upper curve). 



Figure 4. 

The confidence regions for Q rms -ps an d M at 68% CL for Gaussian statistics. 
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